Renormalization scheme for a multi-qubit-network 
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We present a renormalization scheme which simplifies the dynamics of an important class of in- 
teracting multi-qubit systems. We show that a wide class of M + 1 qubit systems can be reduced 
to an equivalent n + 1 qubit system with n > 2, for any M. Our renormalization scheme faith- 
fully reproduces the overall dynamics of the original system including the entanglement properties. 
In addition to its direct application to atom-cavity and nanostructure systems, the formalism of- 
fers insight into a variety of situations ranging from decoherence due to a spin-bath with its own 
internal entanglement, through to energy transfer processes in organic systems such as biological 
photosynthetic units. 



Many-body problems are very difficult, if not impossi- 
ble, to solve exactly. Any simplifications are therefore of 
great potential importance - not only because of prac- 
tical applications but also because of a basic theoretical 
interest. Such simplifications usually result if some exact 
or approximate symmetry can be identified in the under- 
lying Hamiltonian (e.g. Ref. Q). In terms of practical 
applications, most exact results to date have concerned 
systems where all the interacting objects (e.g. particles) 
are indistinguishable (e.g. Ref. [2|). This is understand- 
able, since systems such as many-electron gases have been 
of great experimental interest over the past few decades. 
However given the current levels of activity in the field of 
quantum information processing, there is a clear desire to 
develop such theoretical results for multi-qubit systems. 
In particular, it would be highly desirable to obtain exact 
or approximate results for the generic situation in which 
a collection of qubits interacts with an auxiliary system 
such as a cavity mode or a central spin p|. 

Here we show that a system of M + 1 qubits whose in- 
teractions resemble a spin-star configuration (see Fig. ^) 
can be mapped onto an equivalent n + 1 interacting sys- 
tem with n > 2, preserving the dynamics of the central 
qubit and the quantum correlations of the original sys- 




FIG. 1: Schematic representation of the renormalization 
scheme in which an interacting M + 1 multi-qubit with one 
excitation can be reduced to an equivalent n + 1 qubit system 
with one excitation, where n > 2. For example, the excitation 
can correspond to having a central qubit in its excited state 
(spin up) while the outer qubits are in their ground state 
(spin down). A wide range of interactions are possible be- 
tween outer qubits, including nearest-neighbor, dipole-dipole, 
and pairwise interactions between any pair in the ring. 
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FIG. 2: Schematic diagram of the light-harvesting complex 
LH-I and reaction centre RC structure in purple bacteria. For 
details see Ref. Q . This photosynthetic unit closely resembles 
our model system in Fig. 



tern. This setup can be realized using multi-atom-cavity 
systems, and could also be engineered from a collection 
of quantum dots in an optical cavity. It also mimics 
naturally-occuring photosynthetic complexes (see Fig. |5} 
which are of fundamental importance in nature. More- 
over, we will show that it yields insight into the decoher- 
ence properties of spin-baths possessing internal entan- 
glement. 

The key physical feature which underpins this equiva- 
lence, is a particular non-trivial symmetry in the inter- 
actions experienced by all outer qubits. This symmetry 
allows us to describe the system's dynamics in terms of 
two time-dependent variables, associated with the outer 
qubits and central qubit (or cavity mode) respectively. 
We consider M identical qubits which interact among 
themselves and with a central qubit via the following in- 
teraction Hamiltonian (c.f. Fig. 1): 

M 
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where ^ is the usual raising (lowering) operator for 
j'th outer qubit, or the central qubit with j = C. The 
formalism and results in this paper apply to a wide 
range of possible two-body interactions Jjk, e.g. short- 
range (i.e. nearest-neighbor), long-range (i.e. pairwise, 
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between any pair in the ring) and dipole-dipole. The 
Hamiltonian Hi preserves the number of excitations, i.e 
[H T ,N] = with TV = Y^Li vt a j +< J c (J c- We focus 
hrst on the dynamics within the single-excitation sub- 
space. A basis is given by states in which one qubit 
is excited and the rest are in their ground state, i.e. 
{|«!>|0o>,|©>|0o),...|«m>|0o),|0b)|1o>} with 

\ qj ) = |0 1 ,0 2 ---l J ---0 Aff ) (2) 
|0 fl ) = |0i,0 2 ---0 r --0 M ) 

with j = 1, . . . M. The state of this M-plus-central qubit 
system is given by the unitary evolution associated with 
Hf. 

|*(i)> = |lB>|0 c ) + M*)|0fl)|lo) (3) 

with |l fl ) = J2fLi bj(t)h) ^d E^i \bAt)\ 2 + \bc(t)\ 2 = 
1. We have written the state in this way in order to 
emphasize the collective behaviour of the outer qubits 
- indeed, as a reference to one of the potentially im- 
portant applications of this work, we will frequently use 
the term bath to denote this collection of outer qubits. 
The system's state satisfies the Schrodinger equation 
d\$>{t))/dt = -li?/ 1 which leads to a set of first- 
order coupled differential equations for the complex am- 
plitudes bj(t) and bc(t). Due to a collective symme- 
try in the effective interaction experienced by all outer 
qubits, it is possible to describe the system's dynam- 
ics in terms of two time-dependent variables b c (t) and 
B(t) = Y^jLx bj{t) associated with the central qubit and 
all outer qubits respectively. The nature of this symme- 
try is as follows: Let us assume that each one of the outer 
qubits interacts with the center via an identical coupling 
7jC = 7. Because of this, the strength of the effective 
interaction between a qubit j and the rest of the outer 
qubits is captured by a parameter Aj = J2k^j Jjk> i- c - 
the sum of all the coupling strengths between qubit j 
and any other outer qubit. When Aj is identical for 
all qubits, i.e. Aj = A, as it is for example in the 
case of nearest-neighbour or pairwise interactions, the 
system's dynamics can be described in terms of b c (t) 
and B(t). Note that this does not imply that the pair 
couplings Jjk need to be identical for all possible pairs. 
Instead it implies that the collective interaction expe- 
rienced by one of the outer qubits, due to the rest, 
should be identical for all qubits. Hence the symmetry 
is such that the interaction experienced for each qubit 
(central or outer) is of a collective nature. Given this 
symmetry, b c (t) and B(t) satisfy the following set of dif- 
ferential equations: b c (t) + iAb c (t) + AMj 2 b c (t) = 
and B{t) + iAB(t) = -iM~fb c (t), with solutions of 
the form b c (t) = fo(t)B(0) + g c (t)b c {0) and B(t) = 
f B (t)B(0) + ss(t)MO) where 

f c (t) = -ie- lAt / 2 2 1U (t)/Q 
9c(t) = e- iAt J 2 [(iA/n)u(t)+v(t)] 



f B (t) = e- iAt / 2 [(iA/n)u(t)-v(t)] 

g B (t) = -ie- lAt/2 2M 1U (t)/n (4) 

with u(t) = sin(fh/2), v(t) = cos(Ot/2), and 

Q = n M (rf, A) = ^4Af 7 2 + A 2 (5) 

is a frequency that captures the collective character of the 
interactions in the system. Note that 2\fM^ captures the 
collective features of the interaction between the outer 
and central qubits, while A captures the collective inter- 
action between one outer qubit and the remaining M — 1 
qubits. The relevant expectation values for the dynam- 
ics of the central qubit are determined by the reduced 
density operator p c (t) = Tr outer {|*(t))(^(t)|}. For an 
initial state of the same form as Eq.@ we have 

Pc(t) = \b c (t)\ 2 \0 c )(0c\ + (1 - |Mt)| 2 )|l c )(l c | • (6) 

Hence all the physical properties of the central qubit are 
determined by bc(t) and the relevant expectation values 
can be calculated, e.g. (cr?) = 2\bc(t)\ 2 — 1. 

The system's dynamics are characterized by two ef- 
fective interaction strengths 7 and A. This suggests that 
the entanglement properties should also be describable in 
terms of two contributions [5j : one corresponding to the 
entanglement among all outer qubits, which we call £b, 
and another corresponding to the entanglement between 
the central qubit and the rest, which we call £bc- The 
time-dependent versions of these quantities are given by 

£ B (t) = \-l + \B(t)\ 2 + \b c (t)\ 2 \ (7) 

and 

£ B c(t)=Mb c (t)\ 2 (l-\b c (t)\ 2 ) (8) 

We now take advantage of some known results for 
W— states since these are the states of interest in our 
multi-qubit system (see Eq. [31 . In particular, for any 
partition in the system the entanglement is entirely com- 
posed of pairwise contributions |5( that can be quantified 
by the concurrence 0. The concurrence of the reduced 
state of two qubits in our system has the form Cjk — 
2\bj(t)bk(t)*\ < 1. Hence a measure of the total entangle- 
ment among the outer qubits is E B — k) ^jk, while 
the total entanglement between the central and outer 
qubits is E B c = YljLi Cjc- We now demonstrate that 
£ B (t) and £ B c{t) are lower bounds of E B and E B c, re- 
spectively i.e. £ B (t) < E B , and £ B c(t) < E B c, and 
hence they can be used to quantify the entanglement 
properties of our system. We express the complex ampli- 
tudes bj(t) as bj(t) = \bj(t)\e l0k hence 

M M 

j=i k=i 

M 

i'=l U,k) 
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Using 1 = J2 'Li l&i(*)l 2 + IM*)I 2 , we have 



S B (t) 



CjkCoa(a jk ) 



< Y,C jk \cos(a ]k )\ < J2 

(j,k) (j,k) 
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Knowing that Cjc < 1, we then have that E B c > 
E£i C] c = 4\b c (t)\ 2 (l~\bc(t)\ 2 ) = £ BC {t) which com- 
pletes our proof. 

An important observation from the above relations is 
that the central qubit dynamics (Eq.JfJJl) and the central- 
outer qubit entanglement (Eq. are both completely 
determined by the probability of having the excitation on 
the central qubit, i.e. |&c(i)| 2 , while the intra-bath en- 
tanglement £ B (t) depends on both |6c(t)| 2 and \B (t)\ 2 . 
Therefore |&c(i)| 2 and £ B (t) are the relevant quanti- 
ties to characterize the dynamics of the M + 1-qubit 
system. Henceforth we shall refer to the quantities as 
P M (t) = | b c (t)\ 2 and E M {t) = £ B {t) where the sub- 
script M indicates that these quantities correspond to 
the M + 1 system. 

The above considerations lead to the following dynam- 
ical equivalence, which we shall now discuss: A system 
of M + 1 qubits with interactions forming a spin-star 
configuration and characterized by two parameters 7 and 
A, is dynamically equivalent to a system of n + 1 qubits 
in a similar configuration and characterized by A n and 
j n given that their collective frequencies are identical 
4M7 2 + A 2 = 4n7 2 +A 2 . Since 7 and A are independent 
of each other, we can fix 7 in the equivalent and original 
systems, i.e 7 = j n , such that 



A 2 = 4(M - n) 7 + A 2 



(9) 



The above statements imply that we can find transfor- 
mations between the dynamical quantities of the origi- 
nal and equivalent systems, i.e. Pm (t) = T(P n {t)) and 
EM(t) = F{E n (t)). These transformations are: 



P M (t)-P M {0) = P n (t)-P n (0) 

cum a n 
E M (t)-E M (0) E n (t)-E n (0) 



M 



(10) 



with a n = \B n (0)\ 2 - nP n (0) and /3 n = n{n - l)P n (0) + 
(n - 1)|S„(0)| 2 . Notice that /?„ sets the minimum n m 
the above relations to two, i.e. n m i n = 2. Therefore 
in order to preserve the entanglement properties of the 
original system, the equivalent system should have at least 
two outer qubits. Hence as far as the dynamics of the 
central qubit and the overall entanglement properties are 
concerned, it is possible to map anM+1 qubit system 
onto a system of 2 + 1 interacting qubits (c.f. Fig. 1). 
As an example, consider the case where the excitation 
is initially on the central spin. The time-evolution for 
PM(t)/M and Ej^if) in an original system with M = 10 
and the equivalent system M = 2, are shown in Fig. 2. 
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FIG. 3: (a) Dynamics of central qubit, and (b) intra-bath 
entanglement properties for original and equivalent systems 
when the excitation is initially on the central qubit, i.e. 
6c (0) = 1. Nearest-neighbour interactions have been assumed 
such that A = J and Q, = a/407 + J- We take J = 7. 



They are described by the following relations: 



P M (t) = 
Em (t) = 



1 - 



M 

T 



M 



M(M - 1) 



£ 2 (t) • 



(11) 



We now discuss the effect of dissipation in the system 
using the quantum jump approach[6|. We assume iden- 
tical decay-rates for the outer qubits T, but these can be 
different from the decay-rate k of the central qubit. The 
non-unitary dynamics conditioned on 'no-loss' of excita- 
tion is given by H cond = Hi—iT J2j=i,m a t a j +^c a C • 
The unnormalized state of the system is \^ C ond{t)) — 
exp[— iH con dt]\^(0)) which can be expressed in the same 
form as in Eq.©. Hence, we can again describe the sys- 
tem's dynamics in terms of (normalized versions of) be (t) 
and B{t). In this case, one obtains: 



fc{t) = -ie- xt ' 2 2 1U {t)/n 
9o(t) -- 



-Xt/2 



(6 + 1A) 



fs(t) 

g B (t) 



- p ~Xt/2 



{S + iA) 



fl 



u{t) + v(t) 



u(t) - v(t) 



-ie- xt/2 2M~fu(t)/£l 



(12) 



where X = K+T+iA, S = T— k is the effective dissipation 
coefficient, and u(t) and v(t) are defined as before. The 
collective frequency is now given by 



= M ( 7 , A, S) = ^4M7 2 - (S + i^) 2 



(13) 



With the above equations, one can obtain a similar renor- 
malization scheme to the one discussed earlier. 

The above formalism can be immediately applied to 
model a quantum spin-bath with intra-environmental 
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coupling. The spin-bath model has typically been em- 
ployed in the literature to describe the system-bath inter- 
action at low temperatures. In order to achieve tractable 
approaches, interactions among the bath spins have usu- 
ally been neglected. However, recent works have started 
to explore whether the interaction among bath spins 
might indeed have significant knock-on effects on the dy- 
namics of the central qubit 0, . These studies suggest 
that the intra-bath interactions can suppress decoherence 
of the central spin. It has also been argued that such an 
effect is due to the fact that the intra-bath entanglement 
limits the spin-bath entanglement, and hence the deco- 
herence of the central spin 0. These features can be 
easily understood within our approach, as follows. The 
intra-bath interaction is represented by the effective cou- 
pling A, which affects the dynamics of the central spin 
through the collective effective frequency O. The relation 
between the spin-bath entanglement (£bc(*)) and the 
intra-bath entanglement (£s(t)) can be seen explicitly in 
Eq.(JZj) where it is clear that £s(f) limits |6c(i)| 2 , thereby 
limiting the spin-bath entanglement. Most importantly, 
our approach indicates that such complex many-body 
features can be simulated by a simple system of only 
three interacting qubits. 

Our theoretical formalism can also be used to inves- 
tigate the excitation transfer between a light-harvesting 
LH-I complex and the reaction centre RC in photosyn- 
thetic bacteria (see Fig. 0and Ref. 8]). The LH-I ring 
is made up of 32 donor units - each of these H, in ad- 
dition to the RC, can be treated as a two- level system 
to a good approximation. The prevalent interaction be- 
tween donors is given by an induced dipole-dipole cou- 
pling which can be approximated as Jjk — J/ r jk where 
Tjk is the relative position vector between the outer qubits 
j and k. To a good approximation, all the induced dipole 
moments can be taken as identical and lying perpendic- 
ular to the plane containing the outer qubits, yielding 
A = J^2^. I =2 (^/ r ik) 3 - Hence, as far as the dynamics of 
the center is concerned, this complicated 32-donor ring 
can be accurately represented by two donors. We do not 
pursue this any further here, but leave it as an interesting 
consequence of the present theoretical study. 

We now discuss the extension of these results to larger 
numbers of excitations in the multi-qubit system. Since 
the central qubit behaves as a spin-i particle, only two 
dimensions of the bath's state-space are required to ex- 
pand the pure state in its Schmidt decomposition. For 
any number of excitations N < M, we can express the 
system's state in an analogous way to Eq.JSJ 

\9{t)) = \N B )\0c) + \N-l B )\lc) (14) 

where \Ab) is a time-dependent superposition of (^) 
states, each one having A excitations. We can there- 



fore justifiably claim that the dynamics of a system with 
N < M excitations is analogous to the case of M + 1 — N 
excitations. A particular case is N = M whose solu- 
tions are analogous to the case of single excitation case 
we have discussed. This can be seen in Fig. ^but now 
interpreting the presence of an excitation as a spin-down. 
Therefore, a system with M + 1 qubits and M excitations 
can be mapped on to a 2 + 1 system with 2 excitations. 
The complexity of the dynamics is, however, highly non- 
trivial for M — N > 0. Although we don't currently have 
a full solution, insight into this problem can be gained 
by analyzing the effect of Hj on \^(t)) in detail. In par- 
ticular, we re-write Hi = Vb + Vbc where Vb is the 
interaction between central qubit and the outer qubits, 
and Vbc represents the interaction among outer qubits. 
Vbc induces transitions from \Nb)\0c) to \N — ls)|lc) 
while Vb just produces internal transitions in each \Ab), 
since it preserves the internal number of excitations in the 
outer qubits. Hence we conjecture that for any N < M, 
the dynamics can be described in terms of two collective 
variables: one associated with \Nb)\0c) and the other 
with \N — 1b)|1c)- The open question for future study 
then becomes: What is the minimum number of qubits 
needed to represent such collective properties? 

This paper has presented new equivalence relations in- 
volving multi-qubit systems. In particular, this formal- 
ism provides a new way of simplifying the system-bath 
interactions in open quantum systems. 
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